Kinetics of stimulated polariton scattering in planar microcavities: Evidence for a 
dynamically self-organized optical parametric oscillator 
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We demonstrate for the first time the strong temporal hysteresis effects in the kinetics of 
the pumped and scattered polariton populations in a planar semiconductor microcavity under 
a nano-second-long pulsed resonant (by frequency and angle) excitation above the lower polari- 
ton branch. The hysteresis effects are explained in the model of multi-mode scattering when 
the bistability of the nonlinear pumped polariton is accompanied by the explosive growth of 
the scattered polaritons population. Subsequent self-organization process in the nonlinear po- 
lariton system results in a new — dynamically self-organized — type of optical parametric oscillator. 



PACS numbers: 71.36.+C, 42.65.Pc, 42.55.Sa 



A giant stimulated polariton-polariton scattering is 
one of the most striking features in the optical response 
of planar microcavities (MCs). The scattering was firstly 
observed in GaAs-based MCs with InGaAs quantum 
wells (QWs) in the active layer under a cw excitation at 
wave-vector k p close to the inflection point of lower po- 
lariton (LP) branch WLp(k), when the scattering exhibits 
an unusually low (smaller than 400 W/cm 2 ) thresh- 
old m s, m. Specifically, such excitation results in the 
strong parametric scattering into states positioned ap- 
proximately on WLp(k) with k s = and ki = 2k p , called 
signal and idler, respectively. The effect was theoreti- 
cally described in terms of four- wave mixingorparamet- 
ric scattering 0, [f| [||. Subsequent studies 0, H, 0] have 
shown that the shift of the excitation from the inflection 
point of the LP dispersion is not followed by the corre- 
sponding shift of the stimulated scattering along the LP 
branch characteristic for the four wave mixing. Instead, 
the scattering goes on to the same states with k ~ 
and ~ 2k p . The energy conservation is then fulfilled by 
the shift of the signal and the idler much above the LP 
branch. 

The stability analysis of the single macro-occupied 
pump mode as well as numerical simulations of the po- 
lariton scattering indicate [t| [l(| 11, [12] that such un- 



usual behavior can result from the interplay between two 
instabilities in the resonantly excited MC: bistability of 
the pumped polariton mode intensity with respect to the 
external pump, and its parametric instability with re- 
spect to the decay into multiple scattered polaritons in a 
wide range of k. 

The single-mode optical bistability in MC for pump- 
ing at k p = has been observed and explained theo- 
retically within one-mode optical parametric oscillator 
(OPO) model [3]. The bistability of the scattered sig- 
nal at pumping near inflection point was also found and 



explained within three-mode OPO model [14 1. 

Different possible regimes of the above-threshold OPO 
have been anal yze d theoretically within the three-mode 
approximation [l5l[l6|. However three- mode OPO model 
cannot determine the signal and idler wave vectors k. s j 
which are selected by the parametric process above 
threshold. In the system with a specific dispersion (con- 
taining an inflection point), the interplay between the 
pump single mode bistability and its multi-mode para- 
metric instability can result in a regime where essentially 
multi-mode coupling between the ensemble of lower po- 
laritons plays the decisive role in the formation of the 
signal and idler [ll], [TtJ ■ We will refer to this model as 
a dynamically self- organized (DSO) OPO. 

Dramatic change of LP scattering pattern from 
the figure-of-eight shape corresponding to spontaneous 
regime at lower pump intensity 0, EH to that directed 
along k ~ and ~ 2k p for pump powers above the para- 
metric scattering threshold observed recently [l9| sup- 
ports the DSO OPO model. However, the direct evi- 
dence of multi-mode self-organized nature of the scatter- 
ing demands the study of the LP system dynamics. The 
previous time-resolved studies of the LP dynamics under 
ps-long excitation pulses [IH have revealed a well pro- 
nounced figure-of-eight distribution of the final LP scat- 
tering states in low excitation regime, but only a small 
narrowing of the LP momentum distribution with time 
at higher excitation. 

In this Letter we report the optical studies of kinetics 
of the LP system under ns-long pump pulses and dis- 
cover for the first time the strong hysteresis effects in the 
kinetics of the optical response of the pumped as well 
as of the scattered MC polaritons, thus directly demon- 
strating the behavior, predicted within the DSO OPO 
model IE ESQ- 

The MC structure has been grown by a metal or- 



2 



ganic vapor phase epitaxy. The Bragg reflectors are com- 
posed of 17 (20) repeats of A/4 Alo.13Gao.s7As/AlAs lay- 
ers in the top (bottom) mirrors. The 3/2 A GaAs cavity 
contains six 10-nm thick In 6Gao.g4As/GaAs quantum 
wells (QWs). The Rabi splitting is ft ~ 6 meV. A grad- 
ual variation of an active layer thickness along the sam- 
ple provides a change in the photon mode energy Eq 
and, accordingly, in the detuning A between the exciton 
hojx(k = 0) and photon ftwc(k = 0) mode energies. Ex- 
periments have been carried out on several spots of the 
same sample with A in the range from -1.5 to -2 meV. 

The sample was placed into an optical cryostat with 
controlled temperature. A pulsed Ti-sapphire laser with 
pulse duration of ^ 1 ns, line full width at half maximum 
(FWHM) of - 0.7 meV, and pulse repetition of 5 kHz 
has been used for the excitation of the MC at the angle 
14° relative to the cavity normal. The pump beam has 
been focused onto a spot with diameter of 100 /jm. The 
kinetics of angular distribution of PL signal 7(k, t) from 
the MC has been detected in a wide solid angle around 
the cavity normal by the streak camera with spectral, 
angular, and time resolution of 0.28 meV, 1°, and 70 ps, 
respectively. The transmission signal Ztr (k p , t) has been 
detected with the same streak camera. 

The pump pulse profile Ip(t) is shown in Fig. QJ, 
(dashed curve). The magnitude of Ip(t) determines the 
intensity of the external electric field outside the MC 
\E ext (t, k p )| 2 cx Ip{t). The intensity builds up during the 
first 100 ps and then decreases monotonically (by about 
3 times to t = 1 ns). The pulses with a circular (<r + ) 
polarization and a spectral FWHM of ~ 0.7 meV excite 
the MC with A = -2 meV about 0.5 meV above the LP 
dispersion branch at k p = (k px ,k P y) = (1.96,0) ^m" 1 . 
Figure Q}) shows the recorded kinetics of MC emission 
normal to its plane, I(k = Q,t). The spectra are recorded 
in the <j + polarization. The emission in o~ is about two 
orders of magnitude smaller. 

The MC emission I(k = 0, t) is proportional to 
£ , Qw(k = 0,i)| 2 , the intensity of the k = har- 
monic of the electric field on QW inside MC. Figure [TJd 
clearly shows that 1(0, t) and thus the time dependence of 
\EQw(0,t)\ 2 differs significantly from the exciting pulse 
shape. At low P = 11.5 kW/cm 2 , the signal reaches 
its maximum only slightly later (by ~ 50 ps) than the 
pumping pulse and then decreases quickly, by one order 
of magnitude at t ~ 0.35 ns when the pump intensity is 
still about 60% of the maximum. At P > 15.5 kW/cm 2 , 
the signal behavior changes drastically. After its marked 
decrease (together with the pump) in the range of t = 
0.2 — 0.35 ns the signal starts to grow and reaches the 
second maximum at t ~ 0.85 ns already at the excitation 
pulse fall off. The intensity of this maximum grows by 
more than two orders of magnitude in the range of P 
between 14.9 and 17.2 kW/cm 2 , i.e., threshold-like. 

The recorded time dependences of the MC transmis- 
sion at the pump angle itr(k p ,i) are shown in Fig. [T^, 
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FIG. 1: (Color online) (a) Time dependences of the pump 
pulse Ip(k p ,t) and of the transmission itr(k = k p ,t) at dif- 
ferent excitation densities, (b) Kinetics of the LP emission 
intensity 7(k = 0, t) at different excitation densities. Num- 
bers show the peak pump intensity P in kW/cm 2 . 



(solid curves). / t r(k p ,t) is proportional to the inten- 
sity of the k = k p harmonic of the QW electric field 

.Eq-y^kp), t\ 2 inside MC. Again we see clearly that 
itr(k p ,i) and thus |i?Qw(k p )| 2 do not directly follow 
the exciting field |£ cxt (k p , t)\ 2 cx I P (t). At P < 11.5 
kW/cm 2 7 t r and, thus, |i?Qw(k p )| 2 is a monotonous su- 
per linear function of |.E cxt (k p )| 2 on both the up and 
down going parts of the excitation pulse, the maxima of 
Ztr(i) an d Ip(t) and, hence, those of |i?Qw(k p ), t\ 2 and 
|-Bext(k p , t)\ 2 nearly coincide with each other. However, 
the monotonous dependence I tT (l p) becomes distorted 
with increasing P. Figure shows that I tY starts to 
demonstrate a narrow second peak in the range of nearly 
constant exciting field | £7 CJC t (k p ) | 2 at t ~ 0.2 ns: The 
growth of |£ , Qw(k p )| 2 starts at t — 0.12 ± 0.02 ns, con- 
tinues about 0.1 ns, and gives way to its sharp decrease 
at t = 0.23 ± 0.02 ns. The duration of the increase and 
decrease in |i?Qw(kp)| 2 is close to an available time reso- 
lution of our detecting system of 70 ps. The second peak 
of |I?Qw(kp)| 2 grows quickly with P and shifts slightly 
towards the pulse onset. 

Time dependences Ip{t), I(k — 0, i), and 7tr(k p = 0, t) 
in Fig. [1] can be redrawn as implicit functions I tT (Ip) 
and Isihr) presenting the dependences of the inner field 

£ , Qw(k = k p )| 2 on the exciting field |i? ext (k = k p )| 2 and 
of the k = harmonic of the QW electric field on 

£ , Qw(k = k p )| 2 , respectively, for each P. Figures^ and 
b show the resulting dependences at P = 16.2 kW/cm 2 > 
Pthr- Both the experimentally measured dependences of 
the k = k p electric field inside MC on the external field 
^Qw|k=k p (£ext|k=k p ) and of the k = electric QW field 
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FIG. 2: (Color online) Measured (top panel) and calcu- 
lated (bottom panel) dependences of transmission intensity 
on pump intensity and PL intensity on transmission intensity 
at P = 16.2 kW/cm 2 . Circles mark the characteristic points 
with maxima of transmission (at t = t m ) and k = emission 
(approximately 650 ns later). 



inside MC on that at k = k p i?Qw|k=o(-£'Qw|kp) acquire 
jumps and hysteresis behaviour. 

These experimental results find their qualitative expla- 
nation in the framework of the system of a semi-classical 
Gross-Pitaevskii type equation for QW excitonic polar- 
ization V(k,t) and a Maxwell equation for Eqw(k., t) 
in response to the driving external field far from the 
MC [9, EH [13] • According to this theoretical model, the 
dynamics of the stimulated parametric scattering in the 
planar MCs has a following scenario. Its start is initi- 
ated by a single-mode instability of the pumped mode 
at k = k p , which results in the jump of |i?Qw(k p )| 2 and 
transfers this mode into the region of its strong instability 
with respect to the parametric LP-LP scattering at once 
into a large range of k. That provides an explosive growth 
of LP population in a wide k-space region, mainly around 
k = 0, on one hand, and causes the abrupt decrease in 
the driven mode population, on the other hand. The 
formation of an OPO with a three dominating macrooc- 
cupied modes with signal and idler at k = and 2k p 
occurs due to a dynamical self-organization in the multi 
mode scattering and takes a long - hundreds of ps - time. 

The calculated dynamics of the MC optical response 
within 



the model of Refs. [£j, [ll|, [17| at the excita- 
tion slightly above the hard excitation threshold are dis- 
played in Fig. [2}: and d. The system of Gross-Pitaevskii 
and Maxwell equations has been numerically solved for 
experimental- like time dependence of E oxt (t). As seen 
from Fig. [2fc,d, the model demonstrates the hysteresis in 
dependences of In- on Ip (panel c) and of I(k = 0) on 
Jt r (panel d) similar to the experimental ones. The con- 
sidered model takes into account only coherent scatter- 
ing processes (nondiagonal components of the LP density 
matrix) and neglects important scattering processes like 
LP-phonon or LP-free carriers. Nevertheless, it demon- 
strates the hysteresis in a qualitative agreement with the 




0.5 



' 0.5 



lis 



10 



-0.5 
k (urn -1 ) 



0.5 
Intensity 
(arb. units) 



0.5 
k (Jim" 1 ) 



FIG. 3: (Color online) Time dependences of a fc-distribution 
of LP emission at P = 4.5 (upper panels) and 10.5 (lower 
panels) kW/cm 2 . Left (a,d) and right (c,f) panels display 
the k x dependence at k y = and fc^-dependence at k x = 0, 
respectively. The thick and thin solid lines in the central 
panels (b,e) show, respectively, the pulse and MC emission 
profiles. 



experiment. 

The hysteresis of the dependence of I tI on I p results 
from the blue shift of LP eigenenergy. The increase of 
an overall LP population shifts the driven mode fre- 
quency towards the pump frequency, and results in in- 
creased transmissivity. Even when the pump intensity 
has been fallen down, the system retains considerable 
signal population that keeps it close to the resonance. 
At the same time, the hysteresis of Eqw(]s. s ~ 0) vs. 
itr oc |i?Qw(kp)| 2 reveals a more complicated nature of 
the studied system. Indeed, the model does not pre- 
sume any characteristic time rather than the lifetime of 
cavity polaritons (~3 ps), hence the long time of the 
signal developing (r ~ 10 2 ps) might come from only 
the collective phenomena caused by numerous inter-mode 
scattering processes. Since the evolution to the "three- 
mode" (k = {0, k p , 2k p }) state involves a lot of modes, 
that state appears as an essentially collective formation. 
Moreover, the eventual state differs from stable "three- 
mode" OPO solution even in the case of stationary pump, 
which may be proved by performing the stability analy- 
sis similar to that discussed in 



151 ]. The actual stabil- 



ity of the three-mode pattern is maintained by the pres- 
ence of numerous weak "above-condensate" modes, so 
the whole system occurs to be highly correlated, i.e., it 
demonstrates a new - dynamically self-organized - type 
of OPO. 

To support this scenario of DSO OPO experimentally, 
the LP scattering dynamics in a wide range of k has 
been measured by recording time dependences of an an- 
gle distribution of MC emission. Panels (a,d) in Fig. [3] 
display fc^-distribution at k y = whereas panels (c,f) 
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display k y distribution at k x 
sample with a smaller P t hr r 



=0 for the point on the MC 
3.75 kW/cm 2 for two exci- 
tation densities P = 4.5 and 10.5 kW/cm 2 , i.e., slightly 
and well above Pthr- Ip{t) an d /(k = 0,t) are given 
in panels (b,e). The k— distribution of the emission is 
symmetric in k y direction but shows a well pronounced 
asymmetry in the direction of exciting pulse (k x ). The 
maximum emission in the very beginning of the pulse is 
at kx ^ . 4 (im and then during t ~ 0.1 ns it shifts 
to k = 0. This behavior is well expected. The phonon- 
assisted scattering dominating at low LP densities cannot 
provide the LP relaxation to k = because of comparable 
magnitudes of LP life time and phonon assisted scatter- 
ing time 2(3, 21 1. The effective LP relaxation at higher 
densities appears due to the onset of LP-LP scatterings. 

Figures [3^ and c show that dynamics at P = 4.5 
kW/cm 2 (~ 20 % above Pthr)- The scattering devel- 
oping after the bistable transition at t <~ 0.4 ns occurs 
into a wide range of wavevectors k x between —0.7 and 
+0.3 /urn -1 , k y between +0.6 ^jltcT 1 . With increasing 
time the signal intensity increases about two times and 
reaches its maximum at t ~ 0.6 ns without any marked 
narrowing in its fc-distribution specific for the stimulated 
polariton scattering under cw excitation. 

The marked narrowing of the signal in the fc-space dur- 
ing the excitation pulse duration of ~ 0.8 ns appears at 
higher P. That is illustrated in Figs. [3Jd and d display- 
ing the LP scattering dynamics at P = 10.5 kW/cm 2 ~ 
3 Pthr- The bistable transition of the driven mode at this 
P occurs earlier, at t ~ 0.1 ns, i.e. nearly at the pump 
maximum and results in a highly enhanced LP-LP scat- 
tering. The development of the stimulated signal near 
the band bottom in this case is followed by a monotonous 
narrowing of the A;— distribution of the LP emission both 
in k x and k y directions. Figure 3] displaying the de- 
pendences of a FWHM of fc-distribution shows that the 
dynamical self-organization of the parametric scattering 
takes a long time: the narrowing of the angle distribution 
takes place in the whole time range of the strong scat- 
tering signal up to t ~ 1 ns. The fc-distribution FWHM 
decreases in this time interval from ~ 1 to 0.7 /zm _1 , 
which is still markedly larger than the FWHM in the 
case of cw excitation (< 0.3 ^m _1 ). These experimental 
results clearly prove the DSO OPO model. 

To conclude, the strong hysteresis effects in the kinet- 
ics of the pumped and scattered polariton populations 
have been observed for the first time in a planar semi- 
conductor MC under a nano-second-long pulsed resonant 
excitation slightly above the LP branch. The hystere- 
sis effects are explained in the model of a hard regime 
of the onset of parametric scattering, when the bistabil- 
ity of the nonlinear pumped LP mode is accompanied 
by the explosion-like growth of the scattered LP popula- 
tion and subsequent dynamical self-organization process 
in the open polariton system resulting in dynamically 
self-organized OPO. 
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FIG. 4: (Color online) Recorded kinetics of a FWHM of fc- 
distribution of LP emission Ak x (k y = 0) and Ak y (k x = 0) 
at P = 16.2 kW/cm 2 . 
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